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ABSTRACT 

We perform a global linear perturbative analysis, and investigate the oscillation properties of 
relativistic, non-selfgravitating tori orbiting around neutron stars in the slow rotation limit ap- 
proximation. Extending the work done in Schwarzschild and Kerr backgrounds, we consider 
the axisymmetric oscillations of vertically integrated tori in the Hartle-Thorne spacetime. The 
equilibrium models are constructed by selecting a number of different non-Keplerian distribu- 
tions of specific angular momentum, allowing for disc sizes L ~ 0.5 — 600 gravitational radii. 
Our results, obtained after solving a global eigenvalue problem to compute the axisymmetric 
p-modes, indicate that such oscillation modes could account with most observed lower (i/r,) 
and upper {vtj) high frequency quasi-periodic oscillations for Sco X-l, and for some Z sources 
and Atoll sources with Vl > 500 Hz. However, when i/r, < 500 Hz, p-modes oscillations do 
not account for the linear relation Vtj = Avl + B, B / between the upper and the lower 
high frequency quasi-periodic oscillations that are observed in neutron star low-mass X-ray 
binaries. 
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1 INTRODUCTION 

Quasi-periodic oscillations (QPOs) observed in the X-ray spectra of 
binary systems are transient phenomena associated to non-thermal 
states and state transitions. T hose observed at h igh frequencies, in 
the range ~ 200 - 1200 Hz Jvan der Klisl2005l) . are referred to as 
high frequency QPOs (kHz-QPOs). An important feature of such 
kHz-QPOs is that they usually appear in couples consisting of a 
lower and an upper kHz-QPO. The entire set of kHz-QPOs detected 
in black hole binaries (BHBs), i.e. binaries having a black hole as 
an accretor, contains only seven source^, four of which_| show both 
the upper and the lower kHz-QPO. On the other hand, there are ap- 
proximately 20 neutron stars binaries (NSBs), i.e. binaries having a 
neutron star as the acc retor, that show the kHz-QPO phenomenol- 
ogy (see lMillerl u201Ch for a rece nt review about QPOs in NS Bs), 
including both Z and Atoll typ(Q l lLewin & van der Klisll2006f) . In- 
terestingly, these QPO frequencies correspond to the frequencies 
of orbits a few gravitational radii away from a stellar-mass com- 
pact object, which explains why kHz-QPOs have been considered 
a promising tool to investigate the gravitational forces in the strong 
field regime. It should also be noted, however, that in a recent work 
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1 See lRemillard & McClintockl 1200a) for a review on X-ray properties of 
BHBs. 

2 GRS 1915+105, GRO 1655-40, XTE J1550-564 and H 1743-322. 

3 Based on their timing prope rties and X-ray spectra, NS-LMX Bs are clas- 
sified into Z and Atoll sources (Hasinger & van der Klisl 19891) . 



ISanna et al.l feOld) have shown that the NSB J 170 1-462 provides 
an example of a source that, during the same outburst, presents 
spectral and timing characteristic of both Z and Atoll sources. Be- 
cause the kHz-QPOs detected are remarkably different in the two 
spectral states, and because such differences cannot be attributed 
to changes in the gravitational field of the central compact object, 
the authors conclude that the coherence and rms amplitude of the 
kHz-QPOs cannot be used to deduce the existence of the innermost 
stable circular orbit around a neutron star. 

Although it is still not clear whether there is a unique physi- 
cal mechanism responsible for the generation of kHz-QPOs in both 
BHBs and NSBs, a few observational evidences have emerged over 
the years which may indicate the existence of two distinct mecha- 
nisms. 

• First of all, the kHz-QPO peak separation in NSBs is typi- 
call y within 20% of the neutron star spin frequency /^spm, or half of 
that ( Ivan der Klisl2005l) . In particular, sources with u Bp i n < 400Hz 
have Av ~ u s v j n , while sources with ^ sp in > 400 Hz have 
Av ~ ^spin/2 JMuno et aUl200lh . However, IStrohmaver et al.l 
J2008I) have found that the low-mass NSB 4U0 614+091 does not 
show a strong connection between the kHz-QPO frequency differ- 
ence and the neutron star spin frequency (z/ sp in = 414 Hz, while 
Av ~ 320 Hz). 

• Secondly, while in NSBs the frequencies observed may 
vary by a factor 2 in association to changes of the luminos- 
ity, the frequencies in BHBs are much more stable and vary at 
most by 15%, even when the luminosity changes by a factor 
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3 JRemillard & McClintockl l2006h . This effect is cl early shown 



in Fig . 1 by iBelloni et al.l J2005I) or in Fig. 3 by ITorok et al.l 
( 120060 . where the upper frequencies, uu, are plotted versus the 
lower frequencies, vl ; showing that while kHz-QPOs from a sin- 
gle BHB are represented essentially as a point, kHz-QPOs from a 
single NSB are scattered along straight lines, the so called "Bursa 
lines" l lAbramowicz et al] 20071: lT6r6k et al]|2OO60 . 

• An additional peculiar property of twin kHz-QPOs in the four 
BHBs where they have been simultaneously observed is that such 
frequencies appear in couples obeying the r a tio 3 : 2 to a high de- 
gree of accuracy (see Fig.l by ITorok et al.l ( 120051) ). Over the last 
few years there has been an ongoing debate about the existence 
of the same phenomenology among kHz-QPOs in NSBs. Observa- 
tions, in fact, clearly indicate that Bursa lines vxj = Avl + B of 
kHz-QPOs in NSBs are not compatible with a constant 3 : 2 ratio 
{A — 1.5, B — 0), but it remains controversial whether the peak 
at 1.5 i n the distribution of the observed vy Ivl ratios is physical 
or not (IBelloni et alj|2005l: ITorok et alll2006l : IBelloni et al]|2007l ; 
ITorok et al. 20081: iBoutelier et alj|201oh . 



Several models have been proposed over the years to explain 
the physical mechanism responsible for the generation of kHz- 



QPOs in the X-ray spectra of binary systems ( [Miller et alj 19981: 
|9|; lLamb & Miller 2003; Abra mowicz et all 
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Stella & Vietri 1999 

20031 : iRezzolla et alJ|2003l :IZhang 20 04]: I Arras et alJI 20061) . In par- 
ticular, the model proposed bv lRezzolla et al J ( 120030 for kHz-QPOs 
in BHBs suggests the existence of an oscillating torus around the 
black hole, and identifies the time variability in the X-ray spec- 
tra with inertial-acoustic modes (p-modes) of the relativistic thick 
accretion torus which acts as a resonant cavity for the p-mode ax- 
isymmetric oscillations. The frequencies of the fundamental and 
the first-overtone modes appear approximately in a 3 : 2 ra- 
tio, and within the range of kHz-QPOs depending on the spin of 
the black hole and th e size of the torus. Moreover, the model by 
IRezzolla et al J ( 120030 also accounts for the M _1 scaling of the fre- 
quency, where M is the mass of the black hole. 

The main properties of the axisymmetric p-mode oscilla- 
tions of relativistic, non-selfgravitating tori orbiting black holes 
have been investigated in a series of papers either using a lin- 
ear pe rturbative approach bv lRezzolla et al. (2003): M ontero et al.l 
( 120041) (hereafter Paper I and Paper II, respectively) or through 
non-linear general relati vistic hydrodynamic simulations of both , 
non-selfgravitating tori dZanotti et al.l 1200 3t IZanotti et al.l 120051; 



Montero et alj 120070 and also self gravitating tori JMontero et al.l 

2008. 2010). The linear perturbative approach, computationally 



less demanding, allows for a more detailed investigation of the pa- 
rameter space. Overall, those perturbative analysis confirmed the 
results of non-linear hydrodynamic simulations which revealed that 
the lowest-order eigenfrequencies appear in a sequence of small in- 
tegers 2:3:4..., for a wide range of models. It is also worth noting 
that some of these modes ha ve been related in the analysis per- 
formed bv lBlaes et alj J2OO60 to surface gravity waves, at least in 
the slender torus limit. 

However, it remains unclear whether a kHz-QPO model based 
on axisymmetric inertial-acoustic oscillations is flexible enough to 
explain the rich phenomenology of kHz-QPOs in NSBs, and in par- 
ticular to account for the distribution of the upper and the lower 
kHz-QPOs along the "Bursa lines" observed in NSBs. To address 
this question we perform a global linear perturbative analysis of 
equilibrium thick discs around neutron stars, whose external metric 
is described in the slow limit approximation. The set of perturbed 
relativistic equations is reduced to the solution of an eigenvalue 



problem, following the approach described in Paper I and II. In or- 
der to explore exhaustively the deviations from the 3 : 2 ratio that 
axisymmetric inertial-acoustic oscillations display, we consider a 
much wider parameter space than that studied in Paper I and II, 
particularly in terms of the rotation law and the extension of the 
accretion disc. 

The plan of the paper is as follows: in Section[2]we introduce 
the basic assumptions and equations employed in the definition of 
our general relativistic, vertically integrated tori. In Section [3] we 
derive the perturbation equations, and in Section|4]we list the prop- 
erties of the non-selfgravitating equilibrium models studied. In Sec- 
tion[5] on the other hand, we present the results of the global anal- 
ysis, while Section|6]is devoted to a discussion of the implications 
for explaining kHz-QPOs in NSBs. Finally, Section|7]contains our 
conclusions. 

In the following we will assume a signature { — ,+,+,+} 
for the space-time metric and we will use Greek letters /i, v, A, . . . 
(running from to 3) for four-dimensional space-time tensor com- 
ponents, while Latin letters i,j, k, . . . (running from 1 to 3) will 
be employed for three-dimensional spatial tensor components. We 
also adopt a geometrized system of units by setting c = G — 1. 



2 EQUILIBRIUM TORI IN THE HARTLE-THORNE 
SPACETIME 

Equilibrium thick discs around a rotating neutron star are con- 
structed assuming that their self-gravity can be neglected and that 
the background spacetime t akes the form of the Hartle-Thorne met- 
ric JHartle & ThorneHl968l) . which describes the metric around a 
slowly rotating neutron star. Since we are interested in the region 
of the spacetime around the equatorial plane, we use cylindrical 
coordinates (t, zu, 4>, z), and consider only the zeroth-ord er terms 
in the ratio z/zu < IWilsorj[l972l ; iNovikov & Thorndl 19731) . where 
zu is the cylindrical radial coordinate (see Eq. (1) in Paper II for 
this form of the Kerr line element). Then the Hartle-Thorne metric 
in cylindrical coordinates derives from the Kerr metric, in Boyer- 
Lindquist coordinates, by retaining only those terms which are first 
order in the ratio a/zu, where a is the Kerr parameter, and by re- 
placing the spin of the black hole with the angular momentum of 
the neutron star. In this way we obtain 
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where M is the neutron star mass, and J is t he angular momen- 
tum of the neutron star ( IRezzolla et alJuOOll) . which we assume 
to be constant. Then the line element of the Hartle-Thorne metric 
becomes 



ds 



df 
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zu J 
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(4) 



In order to construct hydrostatic equilibrium models of rotat- 
ing thick discs we solve the continuity equation V a (pu a ) = and 
the conservation of energy-momentum, \7 a T a ^ = 0, where the 
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symbol V refers to a covariant derivative with respect to the met- 
ric I0. Here, T afi = (e + p)u a u^ + pg al3 are the components of 
the stress-energy tensor of a perfect fluid, with u a being the compo- 
nents of the 4- velocity, p the rest-mass density, e the energy density 
and p the pressure. 

It is also useful to introduce an orthonormal tetrad carried by 
the local stationary observer and defined by the one-forms with 
components 



w = zuy/A/Adt 
uj z = dz 



u 4, = s/~A{&4> - u)dt)/m 

u> ra = zu/vAdzu . 



(5) 



In this frame, the components of the four velocity of the fluid are 
denoted by u M and the 3-velocity components are defined as 



it ui a u 



zu, z, 



(6) 



We consider a perfect fluid that follows a polytropic equation of 
state (EOS) p = fcp 7 , where k and 7 = din p/ din p are the poly- 
tropic constant and the adiabatic index, respectively. Next, follow- 
ing Paper I & II, we introduce a vertically integrated pressure 

P(w) = I pdz, (7) 



(8) 



and a vertically integrated rest-mass density 



S(U7) 



pdz, 



where H — H(zn) is the local "thickness" of the torus. We further 
assume that P and E obey an "effective" polytropic EOS 



P = /CE 1 



(9) 



so that K, and r = d In P/d In £ play the role of the polytropic 
constant and of the adiabatic index, respectively. 

After the vertical integration, we enforce the conditions of hy- 
drostatic equilibrium and axisymmetry (i.e. assume dt — = d$) 
and simplify the equ ation of energy-mome ntum conservation to a 
Bernoulli-type form JKozlowski et al.ll 19780 



diP 



di lni 



-ASl 



(10) 



e + p v i-su 

where t = —u^/ut is the specific angular momentum, SI — m^ /u' 
is the angular velocity and 

{u)~ 2 - -(gtt + 2Slg t <p + Sl 2 g 4>4> ) . (11) 

After simple manipulations Eq. d 10b can be rewritten as 



di lni 



i-m 



d l n 



(dtgu + 2Qdig t 4, + ^ dtg^) 



(12) 

By using the metric terms of ((4j» into dlOl l. we derive the following 
force balance equation for a non-selfgravitating disc in the Hartle- 
Thorne spacetime 

1 dP M/zu 2 - Quj(ui + SI) 



E + Pduj 



1-2M) +SIzu 2 {2lj-S})' 



(13) 



where E is the vertically integrated energy density, defined in com- 
plete analogy to ((7} and {8j. 



3 PERTURBATION EQUATIONS 

We next perturb the hydrodynamical equations introducing Eule- 
rian perturbations of the hydrodynamical variables with a harmonic 



Table 1. Main properties of the equilibrium models studied. From left to 
right the columns report: the name of the model, the type of specific angu- 
lar momentum distribution, the constant coefficient £ c , the power-law index 
q (cf. eq. 1211 or its range of values, the minimum radial size for the corre- 
sponding sequence of tori, and the maximum radial size for the sequence of 
tori. The angular momentum of the neutron star is set to be J = 0.1 in all 
of the models. 



Model 


l(m) 


lc 


g 


-L-"min 


J-'max 


Al 


const. 


3.60 


0.0 


0.4 


1.9 


A2 


const. 


3.65 


0.0 


0.5 


5.7 


A3 


const. 


3.70 


0.0 


0.5 


10.0 


A4 


const. 


3.75 


0.0 


0.4 


16.8 


A5 


const. 


3.80 


0.0 


0.8 


29.9 


A6 


const. 


3.85 


0.0 


0.7 


71.4 


Bl 


power-law 


3.0 


[0.1,0.15] 


2.5 


374.3 


B2 


power-law 


3.1 


[0.08,0.15] 


2.1 


348.6 


B3 


power-law 


3.2 


[0.07,0.15] 


1.9 


279.6 


B4 


power-law 


3.3 


[0.05,0.15] 


0.8 


318.4 


B5 


power-law 


3.4 


[0.04,0.15] 


2.5 


363.3 


B6 


power-law 


3.5 


[0.02,0.15] 


1.7 


265.3 


B7 


power-law 


3.6 


[0.006,0.15] 


0.4 


327.1 


B8 


power-law 


3.7 


[0.001,0.15] 


1.9 


276.7 


B9 


power-law 


3.8 


[0.001,0.15] 


2.6 


299.1 


BIO 


power-law 


3.9 


[0.001,0.15] 


2.7 


281.6 


CI 


power-law 


2.59 


0.2 


3.3 


160.8 


C2 


power-law 


2.15 


0.3 


5.3 


548.8 


C4 


power-law 


2.19 


0.3 


7.6 


575.1 


C5 


power-law 


2.29 


0.3 


1.4 


550.9 


C6 


power-law 


2.35 


0.3 


1.9 


545.1 


C7 


power-law 


2.39 


0.3 


0.7 


156.6 


C8 


power-law 


1.79 


0.4 


51.5 


593.9 



time dependence of the type 



} I ~ exp(— icrt) 



(14) 



where 5Q = SP/(E + P) and where we have defined the vertically 
averaged velocity perturbations respectively as 



5V™ = — 
2H 



Sv & dz , SV^ = — / Sv't'dz . (15) 

2# J-H 



We assume that the Eulerian perturbations in the metric functions 
can be neglected, i.e. 5g a t, — (Cowling approximation; Cowling, 
1941). While this condition does not hold in general, it represents 
a very good approximation in the case of non-selfgravitating tori. 

To eliminate the imaginary part from the system of equations 
we introduce the following quantities 



5U = i5V" 



8W = 5V* 



(16) 



and after a bit of straightforward algebra, we derive the following 
set of ordinary differential equations 



o-=8U + a—SQ' + 
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(19) 



where a = 1/(V ) 2 , f = TP/E + P, and the index ' indicates 
the derivative with respect to u>. 

Equations <17M19t are the zo- and (^-components of the 
perturbed relativistic Euler equations and the perturbed continu- 
ity equation, respectively. They can be solved numerically for the 
eigenfrequencies and for the eigenfunctions of p-mode oscilla- 
tions of an oscillating vertically integrated thick disc in the Hartle- 
Thorne spacetime. In practice, we solve the system of equations 
(H7M19I) as an eigenvalue problem using a "shooting" method 
(Press et al., 1986) in which, once the appropriate boundary con- 
ditions are provided, two trial solutions are found, starting from the 
inner and outer edges of the disc respectively, and these are then 
matched at an intermediate point where the Wronskian of the two 
solutions is evaluated. This procedure is iterated until a zero of the 
Wronskian is found, thus providing a value for a and a solution for 
8Q, SU, and SW. The numerical method employed here to solve 
the eigenvalue problem is the same as that discussed in Paper I and 
Paper II, where a more detailed discussion can be found. 



4 EQUILIBRIUM MODELS 

First of all, we need fix the angular momentum of the star which is 
defined as J = /fi*, where fi* = 2n/P+ is the angular velocity of 
the star. In geometrized units we have 



t _ ^ 45 



Mq 

M 



7.1363 x 10" 



(20) 



where J45 is the moment of inertia of the star in units of 10 45 g cm 2 , 
which we take as Z45 = 1.0, while Pi is the period of rotation 
in units of sec. For instanc e, the typical Atoll source 4U 1608- 
52 ( Ivan Straaten et alj|2003r) has a mass M — 1.7Mq, a spin pe- 
riod of P — 1.61ms, and the angular momentum, computed from 
Eq. J20b . is therefore J = 0.15. Even assuming a binary system 
with an accreting neutron star that rotates as fast as the fastest 
known millisecond pulsar PSR J1748-2446ad, namely with P=1.39 
msec, and with a canonical mass M — 1.4M0, would yield to 
J = 0.26. 

In Sections \5\ and [6] below we report results obtained after 
assuming J = 0.1. However, we have also solved the eigenvalue 
problem for the case J = 0.2, without finding any significant dif- 
ference in the results, so that our conclusions remain unchanged. 



Next we define the distribution of the specific angular mo- 
mentum I = l(w) within the disc. We consider tori with distri- 
butions of specific angular momentum that are constant in space, 
i.e. £(zu) = const., and also tori with non-constant distributions 
of the specific angular momenturrQ We note that £(zu) = const, 
is a useful mathematical case which leads to analytic initial data, 
while non-constant distributions of the specific angular momentum 
are a more realistic assumption. In the case of £{vj) = const., the 
value of the specific angular momentum must satisfy the condition 
£ ma < I < l-mb, where £m S and t m \, are the specific angular mo- 
menta of the marginally stable and of t he marginally bound orbit in 
the Hartle-Thorne spacetime (see e.g. lAbramowicz et alj d2003h "). 
On the other hand, in the case of tori with non-constant distribu- 
tions of the specific angular momentum, we consider a power-law 
distribution of the type 



l c w q 



(21) 



where both t c and q are positive constants. The power-law angu- 
lar momentum distributions are chosen such that the position of 
the cusp is always located between the marginally bound and the 
marginally stable orbits. The position of the cusp, as well as the po- 
sition of the maximum rest-mass density ro max in the torus, are ob- 
tained by imposing that the specific angular momentum at these tw o 
points coincides with the Keplerian value JKozlowski et al.lll978T) . 
The inner edge of the torus ro^n is determined by fixing the poten- 
tial gap A Wi„ = Win 



AWm = ln[(-ut)in] - rn[(-w t )cus 






(22) 



On the other hand, the outer edge of the torus roout is defined as 
the position at which P — and it is obtained by integration of the 
hydrostatic balance equation d 1 3 fc . Then, for a given distribution 
of specific angular momentum, sequences of tori having the same 
ro max but different radial extents can be constructed by varying the 
potential gap AWi n - 

In order to investigate how the axisymmetric oscillations de- 
pend on the parameters of the discs, we have constructed sequences 
of models, having different radial extents and different distributions 
of specific angular momentum. The main properties of the various 
models considered are listed in Table[T] Models of class A (hence- 
forth models A) are sequences of equilibrium tori with a constant 
distribution of specific angular momentum, while models B and C 
have a specific angular momentum increasing outwards according 
to Eq. ( 1211 ). Unlike models A and C, which correspond to sequences 
of discs with different radial sizes for a given pair of values of £ c 
and q, models B, refer to different disc sequences which not only 
have different radial sizes but also different values of the power- 
law index q for each of the constant coefficients £ c ; that is for each 
set of models from B 1 to B 10, we have constructed discs with val- 
ues of the power-law index q varying between the minimum and 
maximum values listed in the fourth column of Table Q] at inter- 
vals of Aq = 0.005. The last sequence, C8, correspond to discs 
with a distribution of angular momentum having a power-law in- 



dex q ~ 0.4, close to the Keplerian value q^ a 



0.5. Therefore, 



all these models allow for an extensive investigation of the parame- 
ter space in terms of disc sizes and distributions of specific angular 
momentum, varying from constant to almost Keplerian. 



4 Note that lOian et all 120090 considered an alternative ansatz for the dis- 
tribution of the specific angular momentum. 



kHz-QPOs in NS-LMXBs 5 





1 ' 1 ' 1 


1 1 ' 




\ A6 


- 


2500 




- 




; J»\ q=0.03 




2000 


r \ \ 


- 


£ 1500 




: 


*" 


- Bl \\\ 
q=0.135\\\ 


- 


1000 




- 


500 


i -__c;2 ^^^^^ 


5^^ _ 




C7 

i 


^ S==S **^ ==SS °=*~~=^ 




1,1,1 


1 



20 



40 60 

L A g 



80 



100 



1.54 



1.52 




1.48 



1.46 - 



Figure 1. Eigenfrequencies for the fundamental mode of axisymmetric p 
modes for some representative tori with constant and nonconstant distribu- 
tions of specific angular momentum. Each line corresponds to a sequence 
of tori having the same ro ma x but different radial extents L and the solid 
circles correspond to the values of the Keplerian radial epicyclic frequency 
at OTmax in the Hartle-Thorne metric. 



5 RESULTS OF THE GLOBAL ANALYSIS 

The main properties of axisymmetric oscillations of tori in a Hartle- 
Thorne background are analogous to those found for tori in a 
Schwarzschild and Kerr spacetimes (Paper I and Paper II). Over- 
all, a fundamental mode of oscillation and a sequence of over- 
tones are found (collectively referred to as p-mode oscillations) 
which depend on the position of the rest-mass density maximum, 
on the radial size of the disc, on the distribution of angular mo- 
ment um, while they are rather independent of the equation of 
state jRezzolla et al l200llMontero et alj2004 l Zanotti et al.l2003l : 
IZanotti et al . 2005; Monter o et al.ll2007t) . These properties can be 
summarized as follows: 

• The eigenfrequencies of p-mode oscillations increase as the 
radial size of the disc decreases. 

• The fundamental-mode tends to the values of the radial 
epicyclic frequency at the position of the rest-mass density max- 
imum as the radial size of the tori tends to zero. 

• For any radial extent, the model with the largest fundamental- 
mode eigenfrequency has its rest-mass density maximum located 
at the position at which the epicyclic frequency has a maximum. 

• The ratio between the frequency of the fundamental-mode (/) 
and its first overtone (o\ ) for tori with constant distributions of spe- 
cific angular momentum, £(zu) — const., appear approximately in 
a 2:3 ratio. As the size of the tori tends to zero, the ratio o\ // tends 
tooi//~ 1.52. 

• The ratio between the frequency of the fundamental-mode (/) 
and its first overtone (oi) for tori with nonconstant distributions 
of specific angular momentum, £(zu) — £ c ^ q can deviate signif- 
icantly from oi/f ~ 3/2 for very small discs. As the size of the 
disc increases the Oi/f ratio tends to Oi/f ~ 1.44. 

In FigureQ] we show the eigenfrequencies (in units of Hz and 
scaled for a neutron star mass M = 1.6Mq) for the fundamental 



Figure 2. Ratio oi/f as function of the radial size of the disc for models 
with a constant distribution of specific angular momentum, i.e. models Al 
to A6. 
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Figure 3. Ratio oi/f as function of the radial size of the disc for some 
representative models with a nonconstant distribution of specific angular 
momentum, i.e. models B and C. 



mode corresponding to some representative tori of our sample (i.e., 
A6, B9 with q = 0.03, Bl with q = 0.135, C2, and C7) and re- 
ported as a function of the radial size of the disc expressed in units 
of the gravitational radii r g = GM/c 2 . Each line corresponds to 
a sequence of tori having the same ci7 max but different radial ex- 
tents L and the solid circles correspond to the values of the Keple- 
rian radial epicyclic frequ ency at OT ma x in the Hartle-Thorne metric 
JAbramowicz et alj|2003h . As expected for modes behaving effec- 
tively as sound waves trapped in the disc, the eigenfrequencies de- 
crease like L _1 as the radial extent of the torus increases. Note also 



6 Montero, Zanotti 



1300 
1300 
1100 
1000 
900 
800 

N 

E 700 
600 
500 
400 
300 
200 
100 





1 1 ' 1 ' 1 


1 1 1 1 1 1 1 


1 ' 1 //I / 








// A ' 




B models — 








C models - 







Atoll sources * 


t **i*' ■** 




Z sources 


o ^/ 


Sr^ 9 


. 


Sco X-l 


• vuiw 


•*/ 


- 






' ~~ 


- 


o 


*° of'// 


- 




8° //> 




- 


o * 

■ " A a 


Ay/ 

//'/ 


- 




','/ 




_ 


%* * A //''' 


- 


_ 


■ — /j'*/ 


_ 


- 


- ' AP 


. C 


r X-l 




■" /' 








i XTE ,11807-294 


_ 


■ A? 





3:2 ratio 


- 


,X , i , i 


II 


1,1,1, 



100 200 300 400 500 600 700 800 900 1000 
f (Hz) 




100 200 300 400 500 600 700 800 900 1000 
f (Hz) 



Figure 4. Left panel shows a linear plot of the frequency of the upper kHz QPO versus the frequency of the lower kHz QPO with an asterisk for the Atoll 
sources, circles for the Z sources, solid circles for Sco X-l, squares for Cir X-l and triangles for XTE J1807-294. We also plot, in the relevant range of 
frequencies for kHz QPOs, the first-overtone versus the fundamental mode frequencies for all models listed in Table 1 . Right panel shows the second-overtone 
versus the fundamental mode frequencies instead. 



that, as was shown in Paper I and Paper II for tori orbiting around 
a Schwarzschild and a Kerr black hole, the eigenfrequencies of the 
fundamental mode tend to the values of the radial epicyclic fre- 
quency at ttfmax as the radial dimension of the discs tends to zero. 
As their size diminishes, the role of pressure gradients inside the 
disc becomes negligible and the discs effectively behave as rings of 
particles in circular orbits, oscillating with the epicyclic frequency 
at the maximum rest-mass density point. 

A key feature of the axisymmetric p-modes oscillations of tori 
around black holes is that the eigenfrequencies of the fundamental 
mode and the first overtone appear in an approximately 2:3 har- 
monic sequence, although deviations are possible, in particular for 
the nonconstant specific angular momentum case. This feature is 
also present in the case of sub-Keplerian discs in the Hartle-Thorne 
spacetime. In Figure[2]we show the ratio o\/ f as a function of the 
radial size of the disc for models with a constant distribution of 
specific angular momentum, i.e. models Al to A6. As the size of 
the disc decreases the oi/f ratio increases and tends to a value of 
~ 1.52, independently of the constant distribution of specific angu- 
lar momentum. On the other hand, the behaviour of the oi/f ratio, 
as the disc size decreases, is more complex for nonconstant angu- 
lar momentum disc. In Figure[3]we plot the ratio oi/f as function 
of the radial size of the disc for some representative models with 
a nonconstant distribution of specific angular momentum. These 
models belong to the sequences B and C, and the values of the con- 
stant coefficient l c and of the power-law index q are also shown in 
Figure[3] We observe that for small discs, the o%/f ratio decreases 
as the power-law index q increases for a given value of constant co- 
efficient £ c , i.e. there exists a variation of about 30% in the oi/f ra- 
tio for small discs. In particular, Oi/f has an upper limit of ~ 1.52 
(for discs with an almost constant distribution of specific angular 
momentum), and a lower limit of ~ 1.15 for models with a power- 
law index q close to the Keplerian value qu EV = 0.5. On the other 
hand, for large-size discs the oi/f ratio tends to Oi/f ~ 1.44. 



Table 2. Best fit linear parameters obtained by Belloni et al. (2005, 2007) for 
the Atoll sources, Z sources, Sco X-l and Cir X-l. 



Source 



B 



Atoll sources 0.94 ± 0.02 350 ± 15 

Z sources 0.85 ± 0.01 383 ± 8 

Sco X-l 0.73 ±0.01 469 ±7 

Cir X-l 2.34 ±0.47 104 ± 58 



6 IMPLICATIONS FOR KHZ-QPOS IN NEUTRON STAR 
LOW MASS X-RAY BINARIES 

Based on these properties of th e axisy mmetric p-modes oscillations 
of thick discs, iRezzolla et al.l (2003J) proposed a model of kHz- 
QPOs in BHBs that explains the observed frequencies in terms of 
p-modes oscill ations of a small accretion thic k disc orbiting close to 
the black hole l lSchnittman & Rezzolla 2006). This model accounts 
very well for the A/ -1 scaling of the observed frequencies, for the 
observed variations in the relative strength of the peaks, that are 
interpreted as due to variations in the perturbations that the torus 
is experiencing, and for the fact that twin kHz-QPOs in the four 
BHBs show frequencies obeying the ratio 3 : 2 to a high degree of 
accuracy. As discussed in the Introduction, on the other hand, the 
phenomenology of kHz-QPOs in NSBs present peculiar features 
that distinguish them from those detected in BHBs. In particular, 
the upper and lower kHz-QPOs frequencies uu and vl can vary by 
hundred of Hertz along straight lines vtj = Avl + B, with B ^= 0. 
For convenience, we have listed in Table 2 th e best fit linear param- 
eters obtained by Bel loni et all (I2005, 2007) for the Atoll sources, 
Z sources, Sco X-l and Cir X-l, highlighting the linear correlation 
between Vu and Vl, (errors at la significance level). 

Left and right panels of Figure [4] show the frequency of the 
upper kHz-QPO versus the frequency of the lower kHz-QPO for 
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each class of available NSBs. We have indicated Atoll source^j 
with an asterisk, Z source^j with open circles, Sco X-l with solid 
circles, Cir X-l with squares and XTE J 1807-294 with triangles. 
In addition, for all models for which we have solved the eigen- 
value problem and that are listed in Table 1, we have plotted, in the 
relevant range of frequencies for kHz-QPOs, the first-overtone ver- 
sus the fundamental mode frequencies (left panel) and the second- 
overtone frequency against the fundamental mode frequency (right 
panel). As it is shown in both panels of Figure [4] the portion of 
the plot covered by the observational data intersect only marginally 
with the values obtained from the eigenmode analysis. In particular, 
we show on the left panel, that there are discs (mostly B models) 
with fundamental mode frequency / > 500 Hz and a first-overtone 
frequency which can be in agreement with most observed vl and 
v\j kHz-QPO frequencies for Sco X-l, and for some Z sources 
and Atoll sources. This range of fundamental mode frequencies 
(/ > 500 Hz) indicates that the corresponding suitable models 
would be small in size (L smaller than ~ 50r 9 ), and would have 
a nonconstant distribution of specific angular momentum (B mod- 
els). 

On the contrary, we do not find models lying above the line 
with constant slope 3 : 2, which would be needed to explain 
the observed twin QPOs with vl < 500 Hz. As shown in Fig- 
ure Q] the fundamental mode frequency decreases as the size of 
the disc increases or as the distribution of specific angular mo- 
mentum approaches the Keplerian profile. The oi/f ratio tends to 
oi/f ~ 1.44 as the size of the disc increases or to smaller values 
for small discs with (see Figure |3)- This reflects in the tendency, 
that models show, to concentrate towards the 3 : 2 ratio line as the 
fundamental mode frequency tends to zero (left panel of Figure[4]l. 

The possibility that the observed vl and vjj kHz-QPO fre- 
quencies correspond to the fundamental frequency and to the sec- 
ond overtone of an oscillating torus, encounters similar difficulties 
(see right panel of Figure[4]l. Although the area covered by the com- 
puted p-modes oscillations match some of the observed kHz-QPO 
frequencies, particularly those of Z sources with vl in the range 
150 — 500 Hz, the observations of most of the Atoll sources, of Cir 
X-l and of several Z sources remain unexplained. 

Overall, the properties of axisymmetric p-mode oscillations 
of vertically integrated thick discs are such that, in a plot vjj versus 
vl, the first-overtone and the fundamental mode frequency follow 
a straight line for which A may depart from an exact 3 : 2 ratio 
by 30% (in the case of discs with a nonconstant distribution of the 
specific angular momentum), but for which B « 0. 



7 CONCLUSIONS 

We have performed a detailed analysis of the oscillation properties 
of a thick disc (torus) around a slowly rota ting neutron star. Our 
app roach extends previou s investigations bv lRezzolla et al.l ( 120030 
and lMontero et al.l ( 120041) by considering a much wider parameter 
space and by solving the linear perturbative eigenvalue problem 
in the Hartle-Thorne metric. In particular, the rotation law of the 
torus spans the whole range between a constant distribution of the 
specific angular momentum and an almost Keplerian rotation. We 
have computed a fundamental mode of oscillations and a sequence 



5 4U 1728-34, 4U 0614+09, 4U 1705-44, KS 1731-260, 4U 1735-44, 4U 
1608-52, 4U1636-53, 4U 1820-30, 4U 1915-05, XTE J2123-058. 

6 GX 17+2, GX 5-1, GX 340+0, Cyg X-2, CX 349-2. 



of overtones which can in principle be all excited depending on the 
perturbation acted upon the torus. 

We showed there are discs (B and C models) with fundamental 
mode frequency / > 500 Hz and a first-overtone frequency which 
can be in agreement with most observed vl and vu kHz-QPO fre- 
quencies for Sco X-l, and for some Z sources and Atoll sources. 

However, when these results are used for explaining kHz- 
QPOs in neutron star low mass X-ray binaries with twin QPOs with 
vl < 500 Hz, a major difficulty arises. In fact, unlike kHz-QPOs 
in black hole binaries, the upper and the lower kHz-QPOs in neu- 
tron star binaries obey a linear relation vu = Avl + B, with A 
significantly different from 1.5 (e.g. A « 0.94 for Atoll sources, 
A » 0.85 for Z sources and A w 0.73 for Sco X-l) and B ^ 0. 
On the contrary, the computed axisymmetric p-modes, either in the 
ratio oi/f or 02/ f follow a straight line with 0.8 < A < 1.5 and 
with B « for 01//, and with A > 1.5 and B w for o 2 /f. 

Therefore, with the assumptions adopted in this paper, ax- 
isymmetric p-modes oscillations of a thick disc around a neutron 
star do not provide an explanation for the observed twin QPOs 
in neutron star X-ray binaries with vl < 500 Hz. Neverthe- 
less, additional physics should be taken into account to have a 
better understanding of the differences between the Z and Atoll 
sources. For instance, the thickening of the disc due to radi- 
ation pressure, the interaction of the accreting thick disc with 
the magne tosphere of the neutron star, the presence of a mag- 
netic field (Balbu s & Hawlevll 19911). or non-axisymm etric instabil- 
ities dPapaloizou & Pringlelll984l : lKiuchi et al.ll201ll) may play an 
important role. 

It is known that the luminosities of the Z sources are typi- 
cally close to the Eddington luminosity (L ~ £edd), while the 
luminosities of the Atoll sources appear in a lower and broader 
range (L ~ 0.001— 0.5I/Rnn). Moreover. lHasinger & van der Klisl 
( 119890 first suggested that the differences among the two type of 
sources may be due to differences in the mass accretion rates, i.e. 
low and high accret ion rates show the source as an Atoll or a Z 
one. lLinetai] ( l2009h found that the source XTE J 170 1-462 evolved 
from super-Eddington luminosity to quiescence, displaying an evo- 
lution with features of Cyg-like Z, Sco-like Z and Atoll sources, 
supporting the idea that changes in the ac cretion ra t e are r esponsi- 
ble for this secular evolution. In addition. iLin et a l. (2009) pointed 
out that the Atoll stage may be characterized by a constant inner 
disc radius, while the Z stage may have a luminosity dependent lo- 
cat ion of the inner di sc. Similar results have been recently obtained 
bv lDing et al.l J201 II) . who focused on the interaction between the 
NS magnetosphere and the radiation-pressure-dominated accretion 
disc. As the disc is thickened by radiation pressure, the disc gas 
pressure reduces and the magnetosphere expands, thus making the 
inner disc radius increase. 

Overall, there is a growing evidence that the different source 
stages are due to different disc structures, and that the hot regions of 
the flow where the QPO oscillations are generated could be qualita- 
tively different. Interestingly, the oscillation properties of relativis- 
tic, non-selfgravitating vertically-integrated equilibrium tori orbit- 
ing NSs agree better with the data for the Atoll (low accretion rate) 
stage and the Sco-like stage than for the Z (high accretion rate) 
stage. This seems to indicate that both the NS magnetosphere inter- 
action and the thickening of the disc due to radiation pressure may 
play a crucial role, and will be considered in our future research 
for a more realistic interpretation of kHz-QPO from neutron star 
binaries. 
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